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STORAGE  THEORY  AMD  THE  SUPREMA  OF  CERTAIN 
INFINITELY  DIVISIBLE  PROCESSES 

J.  Michael  Harrison 

Graduate  School  of  Business 
Stamford  University 

1.  Introduction 

Let  X = (X(t),  t > 0}  be  a process  with  stationary,  independent 
increments  (an  infinitely  divisible  process)  and 

-E[X(t ) ] = pt,  where  0 < u < °°  , 

Var[X(t ) ] = o2 t,  where  0 < o2<  °°  , 

for  t > 0.  We  assume  that  X is  continuous  in  probability  with  paths 
that  are  right  continuous  emd  have  left  limits.  We  further  assume  that 
X has  no  negative  jumps,  so  the  Laplace  trains  form  of  X(t)  exists  aind 
has  the  form 

(1)  Ete-0^^]  = e^0^  for  a > 0 and  t > 0. 

The  exponent  function  $(a)  can  be  uniquely  represented  in  the  form 

(2)  ®(a)  = ua  + c2  / x 2(e  001  - 1 + ax)  G(dx)  , 

[o,°°) 

where  G(*  ) is  a probability  distribution  on  [0,»)  satisfying 
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! (3)  / x"1  G(dx)  < « , 

(0,») 


cf.  Gnedenko  (1968),  pp.  323-328.  From  (2)  and  (3)  it  follows  easily  that 
$(• ) is  convex  and  strictly  increasing  with 

(4)  4(0)  = 0,  4'(o)  =n  and  4>"(o)  = o2  . 

Thus  for  each  a > 0 there  exists  a unique  oi(a)  >0  such  that 

(5)  ®(u>(a))  = a (a  > 0)  . 

We  define 


M(t)  = sup(X(s):0  < s < t),  m(t)  = - inf{X(s):0  < s < t) 

for  t > 0-  It  is  known  that  -X(t)/t  -»  4 almost  surely  as  t -*  », 
so  M = lim  M(t)  is  almost  surely  finite.  Let  H(* ) be  the  distribution 
function  (d.f. ) of  M.  We  define  the  first  passage  times  (recall  that 
X has  no  negative  jumps ) 

0(x)  = inf(t  > 0*.X(t)  + x = 0)  , x > 0 


The  following  proposition  was  proved  by  Zolotarev  (1964)  using  analytical 
methods.  An  alternate  proof,  using  combinatorial  methods,  is  given  by 
Takacs  (1967),  p.  88.  Yet  another  (and  I think  simpler)  proof  is  easily 
constructed  from  the  well  known  fact  that 

Z(t)  = e-tD(o)X(t)‘at  , t >0  , 

is  a martingale. 
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Proposition  1:  Efe“°^x^]  = for  x > 0 and  a > 0. 


The  central  purpose  of  this  paper  is  to  prove  by  elementary  means 
that  the  Laplace  transform  of  M is  given  by  wx/$(a).  This  beautiful 
result,  which  greatly  generalizes  the  famous  Pollaczek-Khinchine  formula 
of  queuing  theory,  is  also  due  to  Zolotarev  (1964),  and  again  an  alternate 
combinatorial  proof  has  been  given  by  Takacs  (1967 ) , p.  86.  Our  approach 
relies  heavily  on  the  close  relationship  between  the  maximum  process 


(M(t),  t > 0}  and  the 


process  W = {W(t),  t > 0}  defined  by 


(6)  W(t)  = sup {X (t ) - X(s ):0  < s < t}  = X(t)  + m(t),  t > 


Speaking  loosely,  we  say  that  W is  obtained  from  X by  the  imposition 
of  a reflecting  barrier  at  zero.  The  following  observation  has  been 
made  by  Gani  and  Prabhu  (i960)  and  by  others.  It  follows  immediately 
from  the  fact  that,  for  each  fixed  t > 0,  (X(t)  - X(t-s),  0 < s < t) 
has  the  same  distribution  as  (X(s),  0 < s < t}.  (Each  has  stationary, 
independent  increments,  and  the  one -dimensional  distributions  are  the 


Proposition  2.  For  each  t > 0,  W(t)  and  M(t)  have  the  same  distribution. 


From  this  we  see  that  W(t)  ==>  M as  t -♦  »,  meaning  that  H(*  ) 
is  the  limit  distribution  of  W.  In  section  2 we  study  the  Markov  process 
W,  obtaining  a simple  formula  for  the  Laplace  transform  of  E[W(t)]. 
Zolotarev's  result  is  then  easily  gotten  from  the  fact  that  H(*)  is  the 
unique  stationary  distribution  of  W.  Our  analysis  requires  the  assumption 
that  X(t)  has  finite  variance,  a restriction  not  imposed  in  Zolotarev's 
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original  treatment.  In  section  3 we  shew  how  the  general  result  can  be 
gotten  from  that  for  the  case  of  finite  variance. 

2.  The  Storage  Process 

Generalizing  (6),  we  define  the  storage  process  W with  initial 
state  x > 0 by 

(7)  W(t ) = x + X (t ) + [m(t)  - x]+  , t > 0 . 

For  each  t > 0,  let 

Xt(s)  = X(t  + s)  - X(t)  s > 0 , 

mt(s)  = - inf{Xt(u):0  < u < s),  s > 0 . 

With  W ( * ) defined  by  (7 ) , the  reader  may  verify  that 

W(t  + s)  = W(t)  + Xt(s)  + [mt(s)  - W(t)]+ 

for  t > 0 and  s > 0.  Comparing  this  with  (7),  we  see  that  W is  a 
Markov  process  with  stationary  transition  probabilities  and  state  space 
S = [0,»).  in  the  usual  way,  we  denote  by  probability  distri- 

bution on  the  path  space  of  W corresponding  to  an  initial  state  of  x. 

The  corresponding  expectation  operator  is  denoted  by  Ex(').  We  assume 
that  the  underlying  process  X is  defined  on  some  probability  space 
(0,S,P).  Thus  the  distribution  Px  is  induced  from  P by  the  mapping  (7). 
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We  denote  by  E(*)  the  expectation  operator  corresponding  to 
Let  F^ (t ) = Pf0(x)  < t}  and  observe  that  P{m(t)  > x] 
for  t > 0 and  x > 0.  Thus 


E{[m(t)  - x]  + } = / P{m(t)  - x > y)  dy 
0 


= / P(m(t)  > y)dy  = / F (t)  dy. 
x x y 


From  this  and  Fubini's  Iheorem  we  then  obtain 


(8)  / e_Qrt  E([m(t)  - x]+)  dt  = / / e-cxt  F (t)  dt  dy  . 

0 x 0 y 


Now,  using  Fubini's  Theorem  and  Proposition  1,  we  have 

(9)  J e_at  F (t)  dt  = J c_at  f F (du)  dt 

0 y 0 0 y 


= f /e_at dt  F (du)  = i / 


-au 


0 u 


e F 

a n y 


= i p r e -of®  (y)  i _ l e-“(a)y 
a a 


Combining  (8)  and  (9)  gives 


(10)  / e-at  E([m(t)  - x]+}  dt  = i /"  e^(a)y  dy  = e^(a) 

0 


a 


Now  for  x > 0 and  a > 0 we  define 


cf>  (a)  = / e"^  E [w(t)  ] dt  . 

A A * 


P. 

■ F*lt) 


(du) 


/op(a)  • 


Frcm  (7)  and  (10)  it  follows  immediately  that 
(ll)  <Px(a)  = X^a  ~ + e^a^x/oP(a)  * 

If  X(t)  = Y(t)  - ct,  where  Y is  an  additive  process  and  c is 
a positive  constant,  then  w(t)  can  he  interpreted  as  the  content  of  a dam 
at  time  t.  The  starting  state  x represents  the  initial  content,  Y 
represents  the  input  process  to  the  dam,  and  c represents  the  (constant) 
release  rate.  (See  Gani  and  Prabhu  (i960)  or  Takacs  (1967),  Chapter  6.) 
Thus  formula  (ll)  is  of  considerable  interest  in  itself,  giving  us  the 
Laplace  transform  of  the  expected  content  of  the  dam  as  a function  of  the 
initial  content.  It  has  other  uses  as  well,  however.  From  Proposition  2 
we  know  that  E0[w(t)]  = E[M(t)]  for  all  t > 0,  and  thus 


f (a)  = f e-at  E[M(t)]  dt 
0 


for  a > 0 . 


Since  E[M(t)]  ^*e(m)  as  t ? °»,  it  follows  immediately  that 
afQ(a)  ^E(M)  as  a ^ 0.  Differentiating  (5)  twice  and  using  (1),  we 


obtain 


:>(0)  = 0,  co'  (0)  = l/u  and  -co"(0)  = o ^ [\?  . 


Setting  x = 0 in  (11),  letting  a)§0  and  using  L'Hospital’s  rule,  we 

O 

find  that  cxf^(a)  cr  /2u  as  a\0.  Thus  we  have  the  following 


Proposition  3.  E(M)  = f x H(dx)  = a /2m  . 

S 
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Observe  that  W(fl(x))  =0  if  W(0)  = x.  Also,  X has  the 
strong  Markov  property,  cf.  Hunt  (1956).  Since  0(x)  is  a Markov  time 
for  X,  we  then  have  from  (7)  and  Proposition  2 


Px(W(9(x)  + t)  < y)  = PQ{W(t)  < y]  = P{M(t)  < y) 

i 

for  t > 0,  x > 0 and  y > 0.  Since  ff(x)  is  almost  surely  finite  and  ‘ 

P{M(t)  < y)  -*  H(y)  as  t -*  <»,  it  follows  easily  from  this  that 

(13)  Px(W(t)  < y]  ->H(y)  as  t -* °°  for  all  x > 0 . 

«r- 

*V' 

When  a Markov  process  W has  a limit  distribution  H(*)>  independent  of 
the  initial  state,  it  is  well  known  that  H(*)  is  also  a stationary 
distribution  for  W,  meaning  in  our  case  that 


(l4)  f P (W(t)  < y)  H(dx)  = H(y)  for  t > 0 and  y > 0. 

g x 

The  proof  of  this  general  proposition  for  Markov  processes  is  virtually 
identical  to  that  of  the  corresponding  result  for  Markov  chains,  cf. 
Breiman  (1968),  pp.  13^-135.  (Furthermore,  H(’)  is  the  unique 
stationary  distribution,  but  we  do  not  actually  need  this  fact. ) From 
(ll+ ) and  Proposition  3 it  follows  that 

/ E [W(t)]  H(dx)  = <t2/2u  for  all  t > 0. 

S X 

Using  this,  the  definition  of  <px(a),  and  Fubini's  Theorem,  we  then  have 
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(15)  / cp  (a)  H(dx)  = / e-^  f E [W(t) ] H(dx)  dt 

s x 0 s 


= / e_Qt  (ct2/2h)  dt  = a2/2uQ 


But,  directly  from  (11)  and  Proposition  3, 


/ <f>(a)  H(dx)  = a /2\±a  - u/cr  + [/  e' 


H(dx)  ]/ou>(a) 


Combining  (15)  and  (l6),  we  have 


/ e~a^a)x  H(dx)  = pa>(o)/a  . 


Replacing  a by  ®(a),  and  hence  oj(a)  by  a,  gives  the  central  result. 


Proposition  U.  E(e  uw)  = / e”  H(dx)  = pa/$(a)  for  a > 0 . 

S 


3.  The  Case  of  Infinite  Variance 

Suppose  that  X is  as  described  in  section  1 except  that 
Var[X(t)]  = «.  The  Laplace  transform  of  X(t)  still  exists  and  has 
the  form  (l),  but  the  exponent  function  <t>(&)  cannot  be  represented  in 
the  form  (2).  For  each  k = 1,  2,  ...  let  X^  = (Xk(t),  t > 0]  be 
obtained  from  X by  truncating  all  jumps  of  size  greater  than  k 
(replacing  each  by  a jump  of  size  k).  Then  VarfX^t)]  < ».  We  denote 
by  ^(a)  and  the  negative  of  the  mean,  the  exponent  function, 
and  the  supremum  respectively  of  X^.  For  each  t > 0,  X^(t)  ^X(t)  a.s. 
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as  k f oo.  Thus  pk^|  p,  ®k(a)  V 4>(a)  for  all  a > 0,  and  f- M 
a.s.  as  k Combining  these  facts  with  Proposition  4 and  monotone 

convergence,  we  have 


I 


r 

i 

* 


) = uka/4>k(a) 


->  p.a/$(a) 


E(e 


-aM 


) 


as  k ->».  This  same  truncation  argument  can  be  used  to  show  that  formula 
(ll)  also  remains  valid  when  VarfX(t)]  = «. 
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Abstract 

\\  Let  X - {X(t),  t Z>  0}  be  a process  v>ith  stationary,  independent 
increments  and  no  negative  jumps.  Let  W = {W(t),  t-  > 0)  be  this  same 
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process  modified  by  a reflecting  barrier  at  zero  (a  storage  process). 

Assuming  that  -E[X(t)]  = ut  > 0,  let  'M  = sup{X(t):t  > 0.1,  and  denote 

■ i 

by  i>(a)y  the  exponent  function  of  X.^  A simple  formula  is  derived  for 
the  Laplace  transform  of  E[W(tt],  t 0,  as  a function  of  W(0).  losing 
the  fact  that  the  distribution  of  M is  the  unique  stationary  distribution 
of  the  Markov  process  W,  this  yields  an  elementary  proof  that  the  Laplabe 
transform  of  M is  u«/i>(u''.  If  VarfX't)]  = a t < »,  it  follows  that  \ 
E(M’  = o2/^!.  These  surprisingly  simple  formulas  were  originally  obtained 
by  Zolotarev  using  analytical  methods. 
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